Spin Amplification for Magnetic Sensors Employing Crystal Defects 
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Recently there have been several theoretical and experimental studies of the prospects for magnetic field 
sensors based on crystal defects, especially nitrogen vacancy (NV) centres in diamond. Such systems could 
potentially be incorporated into an AFM-like apparatus in order to map the magnetic properties of a surface at 
the single spin level. In this Letter we propose an augmented sensor consisting of an NV centre for readout and 
an 'amplifier' spin system that directly senses the local magnetic field. Our calculations show that this hybrid 
structure has the potential to detect magnetic moments with a sensitivity and spatial resolution far beyond that 
of a simple NV centre, and indeed this may be the physical limit for sensors of this class. 



A key goal for future sensor technologies is the detection 
of weak magnetic fields at molecular length scales. There are 
numerous potential applications in materials science, medical 
science and biology. An ambitious goal would be to detect the 
field due to single nuclear spins, in order to gain direct infor- 
mation about the structure of a molecular complex deposited 
on the surface. This would require unprecedented combina- 
tion of sensitivity and spatial resolution. 

Among the most sensitive magnetic field sensing devices 
are Hall sensors SQUID sensors [2], force sensors [3], 
and potentially NV centres in diamond [4- 8 1. In this Letter we 
propose a significant improvement of the NV centre sensor by 
coupling it to an amplifying spin system (see Fig.[T]i. We show 
that the improvement in magnetic moment sensitivity can be 
expected to be of about three orders of magnitude; indeed the 
augmented NV centre may have the highest in-principle per- 
formance of any such device. 




FIG. 1. (Color online) (a) The sensor structures: an NV centre (red 
sphere) is embedded in the middle of a diamond nano-crystal that 
is attached to an AFM-tip (yellow cone). The magnetic field gen- 
erated by local spins (green spheres) affects the dynamic of the NV 
centre. The strength of the magnetic field can be inferred by mea- 
suring the NV centre through optical means and manipulating it with 
microwaves, (b) Amplified NV centre sensor: the surface of the nan- 
odiamond is decorated with an another spin system (blue sphere) that 
couples to the NV centre inside the diamond. This additional spin 
system has an amplifying effect and enormously increases the mag- 
netic field sensitivity of the sensor and its spatial resolution. 

We begin with a general discussion of measuring an un- 
known magnetic field using two energy levels of a probe spin 



s system {2s G N), the underlying principle of the NV cen- 
tre sensor Let the probe spin be governed by the Zeeman 
Hamiltonian H = — /ipi.b(Bo +B)S^, where /ipib is the mag- 
netic moment of the probe spin, is the usual spin operator, 
Bo is a known external field applied in the z-direction, and B 
corresponds to the magnetic field (also in the z-direction) that 
we wish to estimate. Consider two eigenstates |0) and |1) of 
H whose z-projections differ by an integer m. The field es- 
timation then proceeds as follows: we start by creating the 
superposition |+) [here |±) — 1 /^/2{\Q) ± 1 1))], for example 
by using a suitable microwave or radio-frequency pulse se- 
quence. This state evolves in time to |'/(f)) = l/v^(|0) + 
exp(/m/Zprb(Bo and therefore by successive mea- 

surements of the spin, we can infer the strength of the mag- 
netic field. In any real-world experiment the state will 
also be affected by decoherence processes, typically dephas- 
ing is predominant |9|. In this case, the coherence between |0) 
and |1) decays as exp(— 7(f)) for a positive non-decreasing 
function 7(f). The time evolution of the two level system's 
density matrix p(r) is then fully described by 

P{t)-^\ , ,„Mp,-b(fio+fi) ^ ^ I • (1) 



The uncertainty of estimating B is limited by the quantum 
Cramer-Rao bound (CRB), i.e. even for an ideal system and 
perfect measurements the uncertainty in 5B cannot become 
smaller than bcK- The CRB is given by the inverse square 
root of the quantum Fisher information F with respect to B 
ifTOl . yielding the following inequality 

8B>bcK^A^ = ^zvT^. , (2) 



/p mTljUprbl 

where T is the duration for which the spin has experienced the 
magnetic field. Repeating the measurement many times gives 
a further reduction in the uncertainty of the unknown field. 
Assuming the preparation of the |+) state takes the time tp, 
the lower bound for the uncertainty 5B after = T /{x + tp) 
repetitions within a total time T is given by 



bcR 
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As an example, we set 7(f) = j^. Minimizing the CRB with 
respect to T yields the optimal sensing time for each run of the 
protocol: 

- ^ (t2 - 2t„ + ^jTi+4tp{3T2+tp)^ . (4) 

It is easily seen that t* « in the limit of weak dephasing 
{T2 ^ tp), in contrast to T* w T2 when the dephasing time is 
comparatively short, T2 ^ tp. 

Since a magnetic field smaller than bcR cannot be resolved, 
the CRB provides a natural limit of the sensitivity of a mag- 
netic field sensor Interestingly, this limit does not depend on 
the actual value of B. However, a naive intuition is that the 
proximity of the probe spin and the sample will be an impor- 
tant characteristic of a practicable sensor, since it determines 
the strength of the measured magnetic field as well as the sen- 
sor's spatial resolution. To quantify this intuition, one can ask 
the question: how long does it take until the sensor spin de- 
tects the magnetic moment associated with an electron or a 
nuclear spin? In the following, we introduce a magnetic mo- 
ment sensitivity (in units of T/y/Uz) as a suitable figure of 
merit for this question: the magnetic moment sensitivity 5 of a 
sensor is given by the number of proton magnetons that can be 
resolved within the time window \/T for a distance r between 
the probe spin and the sample. The dipole field originating 
from a proton with magnetic moment jXp at the position of 
the probe spin is of order b — ^ ^ (where jUq is the vacuum 
permeability), so we obtain 

^^min,(fecR)^ 
b 

For example, once again assuming 7(f) — jr and employing 
Eq. Q gives the following magnetic moment sensitivity 

i — , , — . (o) 

We shall now briefly discuss the example of field sensing 
with an NV centre before describing the inherent limitations 
of this approach and presenting our proposal for improving 
the characteristics of this class of device. An NV centre in 
diamond possesses spin 1 with the three levels |0) and | ± 1). 
The levels | ± 1) are degenerate in the absence of an exter- 
nal field and shifted by a zero-field-splitting (ZFS) of about 
2.87 GHz. For the purpose of sensing the magnetic field, vari- 
ous two level manifolds can be used, for example { 1 0) , | ± 1 ) } 
or{l/V2(|l) + |-l)),(l/V2(|l)-|-l))}0. 

Figure |2] shows the magnetic moment sensitivity of various 
classes of sensors. Clearly, NV centres hold the promise of 
a high field sensitivity combined with a fairly small probe to 
sample separation. However, in order to further improve the 
spatial resolution and obtain an even better magnetic moment 
sensitivity (according to the definition in Eq. the NV cen- 
tre must be brought even closer to the sample. Due to the 
cubic dependence of 5 on r even a modest reduction in the 




FIG. 2. (Color online) Performance of various state-of-the-art and 
proposed magnetic field sensors (figure adapted from Ref. |11|). 
The plot shows magnetic sensitivity (per \/Hz) (horizontal) versus 
a typical tip-sample separation. The data points are reported sensi- 
tivities from experiments: force sensors I12I414I . Hall sensors |15|, 
SQUID sensors |2, 16, 17 1, and NV centre sensors |4|. The diagonal 
lines sketch the boundaries for the magnetic moment sensitivity (see 
Eq. (|5]l) required to detect 1, lO"*, 10*. or 10^ protons (solid black) 
and a Bohr magneton (blue dash dotted) within one second. The po- 
tential of our proposed pre-amplified NV centre is illustrated by the 
dark green ellipse. 



separation leads to a large gain, and for this reason NV cen- 
tres have been embedded in smaller and smaller nano-crystals 
im |6| . However, the size of the diamond crystal surround- 
ing the NV centre cannot become arbitrarily small without 
severely affecting the remarkable coherence time (and thus 
the sensitivity) of the NV centre. In '^C isotopically enriched 
bulk diamond the room temperature coherence time can be as 
long as 1 .8 ms [18] , increasing up to 2.4 ms with dynamic de- 
coupling |fT9l . On the other hand, the best reported coherence 
time for a nanodiamond with a diameter of 7 nm is reduced 
to 1.4 |i.s ||20) . Smaller crystals with a diameter of only 5 nm 
have been studied in Ref. 11211 . however the properties of en- 
closed NV centres are even further degraded in this case. A 
miniaturisation of the crystal hence leads to a reduced field 
sensitivity in exchange for a smaller separation between the 
probe and the sample. 

We propose to overcome this trade-off situation by bringing 
the NV centre "effectively" closer to the sample without re- 
ducing the size of the nanocrystal. This is achieved by attach- 
ing an additional spin system on the surface of a nanocrystal 
which relays the sample magnetic field to the NV centre (see 
Fig.[T]l. In the simplest implementation, a single electron spin 
serves as the amplifier. Of course, further improvements will 
be possible by using higher spin systems such as N@C6o ll22ll 
or a molecular magnet (TT, '24]. It may even be possible to 
make use of more sensitive entangled states with suitably en- 
gineered molecular systems Il25ll26ll . 

For simplicity, we consider an electron spin as the ampli- 
fier in the following. The NV centre and the amplifier spin 
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are dipolar coupled, and we assume the vector connecting the 
two spins is aligned with the z-axis. Both spins experience an 
(known) homogeneous external magnetic field Bq = (0,0,Bo) 
and a small magnetic field whose z-component B at the po- 
sition of the amplifier is to be measured. In comparison to 
the amplifier spin the NV centre experiences a weaker sample 
field cB where the factor c < 1 depends on the relative separa- 
tions. The Hamiltonian thus reads: 



H — -^nv(Bo + cB)5nv,z+I>nv'Snv,z~Ma(5o+5)'5a,z 

+ d{2Sj^y,SA^z~S!^y.xSA.x~S!^V.ySA,y) , (7) 



where d = ^ ^^J"*^' is the dipolar coupling constant; Ynv and 
Ya are the gyromagnetic ratios of the spins, A is the distance 
between the NV centre and the amplifier spin, and D^y de- 
notes the NV centre's ZFS constant. Our protocol assumes 
that the influence of cB on the NV centre is much smaller than 
its coupling to the amplifier spin; this will certainly be justi- 
fied if the amplifier is an electron spin and the sample field is 
due to nuclear spins. In addition we assume that flip-flops be- 
tween the spins are heavily suppressed and can be neglected, 
for example because Dnv represents the largest energy in the 
system, yielding the following effective Hamiltonian 

H W — ^NvBo^NV^z +£'nv5nv^j; 

^ IJ-AiBQ+B)SA,z+^dSm.zSA,z : (8) 
whose level structure is schematically depicted in Fig. [3] 
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FIG. 3. (Color online) (a) Level structure of an NV centre. The 
ground state is a spin triplet with zero field splitting and the excited 
state has a manifold of levels. Following excitation with green light 
the dominant decay is spin preserving. However, spin is not always 
conserved due to an additional decay path to |0) via a metastable 
state, so that continuous excitation eventually results in a polarised 
state, (b) Eigenspectrum of H for an amplifier spin j system. 



The field sensing protocol now proceeds as follows. Ini- 
tially, the system starts in a completely mixed spin state, gig, 
so that it needs to be polarised with laser and microwave 
pulses. As a first step we polarise the NV centre by illumi- 
nating it with green light, after about 1 |j,s the NV centre will 
have relaxed to the state |0) with high probability |5J, leaving 
the combined system in the state |0)(0| (i) ^l2- Next, we ap- 
ply a selective microwave ;r-pulse on the transition |0) | I) ^ 
1 1)1 \.). For a distance A — 10 nm between the NV centre 



and the amplifier, this transition is split from the neighbouring 

|0) 1 1) ^ 1 1) 1 1) transition by 2d w 0.65 Mrad/s 

(i.e. 104 kHz). A highly selective ;r-pulse is beneficial for the 
presently discussed protocol. For this reason we now estab- 
lish the conditions under which imperfections in the pulse se- 
lectivity can be considered negligible. Assuming Lorentzian 
line broadening of both transitions with a FWHM of 2/72_nv 
(72,NV is the coherence time of the NV centre), a detailed anal- 
ysis [271 shows that a linewidth overlap of up to 10% meets 
this requirement and translates into r2 nv > 10 |i.s. The prod- 
uct of pulse duration and frequency passband T x Ao) is typ- 
ically between 2 and 5 |28|. Taking T ~ 4/Aa) and allowing 
the passband to overlap with the area of unwanted transition's 
Lorentzian by at most 5% (allowing us to neglect imperfect 
pulse selectivity), we thus obtain Xy, = 2(2^-cot(0.05;r)/7-,,Nv) 
for the duration of the desired selective ;r-pulse (reducing to 
Tji^ \/d in the limit of long r2,Nv)- Finally, the state of the 
NV centre is measured by exciting it once more with green 
light and detecting the spin-dependent photoluminescence. 
Importantly, this puts the NV centre back into the |0) state, 
leaving the system in one of the well-defined states |0) | \) or 
|0)|t>. 

Without loss of generality we assume the system is in 
the state |0)| J,), ready for the magnetic field estimation 
through several repetitions of the following sensing cycle: A 
microwave | -pulse creates the superposition l/\/2(|0|) + 
|0 t))- It is safe to assume that the amplifier spin can be ro- 
tated fast and with high fidelity as long as the unknown field 
coiTesponds to the smallest energy scale of the system, as can 
be accomplished with a modest external field. After a sens- 
ing time T, the amplifier spin acquires a relative phase pro- 
portional to B. The phase is first mapped onto a population 
difference between |0 |) and |0 t) with another ^-pulse, and 
then entangled with the NV centre spin state with a selective 
;r -pulse in the same way as in the initialisation process. The 
NV centre spin state is now read out using spin-dependent 
photoluminescence detection (also initialising it for the next 
sensing cycle). 

Having described the protocol, we now benchmark the sen- 
sitivity improvement of the amplified system over a conven- 
tional single NV centre sensor Obviously, the performance 
of both sensors depends on the decoherence model, which we 
assume to be fully characterised by y{t). Different forms of 
7(f) corresponding to different predominant dephasing mech- 
anism occur, however, typically 7(f) can be written as {tjT^Y 
for 1,2 or 3 |5 , 20]. 

For our comparison, we consider the levels |0) and 1 1) of an 
NV centre in the middle of a nanocrystal with radius lOnm 
and a control and measurement duration of fp.NV ~ 1 M^s per 
cycle. In contrast, using the electron spin on the surface of 
the augmented sensor entails an additional control overhead 
of about In, i-C- fp.A ~ fp.NV + '^n- Figure |4] shows the ratio 
'5nv/5'a as a function of the amplifier coherence time 72 a 
for several illustrative values of 72.nv- The red regions of 
the inset show the parameter space for which our protocol 



is not necessarily advantageous, either due to excessive line- 
broadening preventing a highly selective ;r-pulse on the tran- 
sition |0)| 4-) i), or because 72^nv is so much larger than 
T2,A that it more than compensates the benefit of the amplifier 
spin. Figure |4] impressively demonstrates that the magnetic 
moment sensitivity can be enhanced by up to three orders of 
magnitude for r2,NV ~ 72,a with realistic parameters. Fur- 
ther improvements may be feasible by substituting the elec- 
tron spin amplifier with a high spin system or molecular mag- 
net. 




FIG. 4. (Color online) Ratio Sj^y/Sa as a function of the amplifier 
spin coherence time T2^a for ?2.NV = 100 \is (black), 30 [is (red), and 
15 VIS (blue). This ratio gives the factor by which the amplified sen- 
sor outperforms a conventional NV centre sensor. Here y(?) =t/T2, 
but we note that the curves for = {t/T2)^"' look almost identi- 
cal. Other parameters: r^t = 1 nm, /"nv = 11 nm, A = 10 nm, fy,.NV = 
1 vis,fp A = 1 Vis + Tji. Inset: excessive line-broadening or a signifi- 
cant mismatch in the coherence times prevent a sensitivity enhance- 
ment in the red region (see text). 

We now turn to another important consideration: Given a 
sensor that can resolve magnetic fields at the level of the single 
Bohr or even nuclear magneton, its spatial resolution is vitally 
important for many applications such as mapping the con- 
stituent molecules and nuclei of complex organic molecules. 
Unsurprisingly, it is highly beneficial to be able to move very 
close to the sample in order to reliably distinguish individual 
magnetic moments. Figure [5] shows the z-component (i.e. the 
measured component) of the magnetic field for a regular ar- 
ray of protons at two different heights above the surface. If 
the probe sample separation is of order the distance between 
individual dipoles, the nuclear spins can be unambiguously 
resolved (given the required field sensitivity). On the other 
hand, the field is almost indistinguishable from that of a sin- 
gle stronger dipole if the sensor is scanning at z = 10 nm, even 
with the ability to sense small magnetic fields of order 1 nT. 

In summary, we have proposed a practical way of enhanc- 
ing N V centre based magnetic field sensor with a pre-amplifier 
system. Such an augmented sensor possesses a magnetic mo- 
ment sensitivity that is capable of resolving individual pro- 
tons, both in terms of their field strength as well as spatially 
on an atomic scale. The sensor consists of two essential parts: 
the NV centre for optical readout coupled to an entirely spin- 
based system for the magnetic field detection in close proxim- 
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FIG. 5. (Color online) The z-component of the magnetic field of a 
regular array of protons calculated in two planes z = 1 nm and z = 
10 nm. The 36 dipoles with the magnetic moment of a proton are 
placed in the z = plane, pointing into the z-direction as indicated 
by the red arrows. The field range of the contour plots spans 43.2 nT 
in the z = 10 nm plane and 2.61 \iT for z = 1 nm. Individual dipoles 
are easily resolved at z = 1 nm, but at a distance of 10 nm their field 
looks similar to that of a single stronger magnetic moment. Spatially 
resolving individual spins at this distance requires a sensitivity that 
goes far beyond the amount of detail reflected in this plot. 



ity to the sample. Despite this partitioning of tasks, the fun- 
damental complexity of this improved sensor design is similar 
to that of existing proposals and implementations. 
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